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FOREWORD 


This  investigation  was  performed  for  the  Commander  and  Director, 
U.S.  Army  Construction  Engineering  Research  Laboratory  (CERL)  under 
Project  4A761101A91D,  In-House  Laboratory  Independent  Research  Program. 

This  investigation  was  performed  by  the  Engineering  and  Materials 
Division  (EM),  CERL,  under  the  direction  of  Dr.  J.  D.  Prendergast,  Prin- 
cipal Investigator,  with  technical  advice  and  support  from  Dr.  W.  E. 
Fisher,  Chief  of  Engineering  Team  of  EM;  Prof.  A.  H.  S.  Ang,  University 
of  Illinois,  Urbana;  and  Dr.  L.  R.  Shaffer,  Technical  Director,  CERL. 

Dr.  G.  R.  Williamson  is  Chief  of  EM  and  COL  J.  E.  Hays  is  Commander 
and  Director  of  CERL. 
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PROBABILISTIC  CONCEPT  FOR  GRAVITY  DAM  ANALYSIS 


1 INTRODUCTION 


Problem 

The  Dam  Safety  Act,  Public  Law  92-367,  directed  the  Secretary  of 
the  Army,  through  the  Corps  of  Engineers,  to  initiate  a program  of 
safety  inspections  of  dams  throughout  the  United  States.  As  a result  of 

these  inspections,  it  was  learned  that  there  are  nearly  20,000  poten- 

tially hazardous  dams  in  the  United  States  which  would  seriously 
endanger  life  and  cause  extensive  damage  to  property  if  they  failed. 

The  recent  failure  of  the  Teton  Dam  and  the  near  failure  of  the  Van 

Norman  Dam  during  the  1971  San  Fernando  earthquake  are  but  two  instances 

during  the  past  several  years  which  have  highlighted  the  need  to  develop 
more  effective  dam  safety  evaluation  techniques  to  yield  the  best  assur- 
ance of  dam  safety  possible  within  the  limitation  of  current  knowledge 
available  in  the  scientific  and  engineering  communities. 

The  Corps  of  Engineers  has  full  responsibility  for  the  safety  of 
nearly  550  dams,  i.e.,  dams  which  the  Corps  plans,  designs,  constructs, 
and  operates. * Most  of  these  dams  are  either  earth  or  gravity  type  as 
shown  in  Table  1. 

In  the  conventional  design  process  for  dams,  many  uncertainties  are 
associated  with  the  forces  e.g.,  external  water  pressure,  uplift,  earth 
pressure,  earthquake  loading,  etc.,  and  the  sliding  and  overturning  re- 
sistance capacities  which  must  be  considered  in  the  design  of  dam  struc- 
tures. Professional  engineers  are  aware  of  these  uncertainties;  how- 
ever, little  attempt  has  been  made  to  quantify  them.  Consequently, 
conventional  designers  rely  heavily  on  the  subjective  choice  of  the 
design  parameters  supplemented  with  a conservative  safety  factor  se- 
lected on  the  basis  of  experience  and  judgment.  Certain  shortcomings 
may  be  observed  with  the  current  design  practice,  namely: 

1.  The  safety  of  the  design  with  respect  to  various  modes  of  fail- 
ure such  as  sliding,  overturning,  and  overtopping  is  unknown,  and  there 
is  not  a quantitative  and  consistent  basis  for  comparing  the  safety  of 
various  designs. 


Dam  Safety  Program  in  the  Corps  of  Engineers,  Report  to  the  Federal 
Coordinating  Council  for  Science,  Engineering  and  Technology  (U.S. 
Army  Corps  of  Engineers,  Office  of  the  Chief  of  Engineers,  1977). 


Table  1 


Dams  for  Which  the  Corps  Has  Full  Responsibility 
[From  Dam  Safety  Program  in  the  Corps  of  Engineers,  Report  to  the  Federal 
Coordinating  Council  for  Science,  Engineering  and  Technology  (U.S.  Army 
Corps  of  Engineers,  Office  of  the  Chief  of  Engineers,  1977)] 


TYPE 

AGE 

DAM  TYPE 

NUMBER 

CONSTRUCTION  DATE 

NUMBER 

Earth 

315 

Before  1900 

17 

Rockfi 1 1 

19 

1900-1909 

8 

Gravity 

141 

1910-1919 

15 

Arch 

3 

1920-1929 

12 

Other 

60 

1930-1939 

70 

1940-1949 

64 

1950-1959 

97 

1960-1969 

156 

_ _ _ 

1970-To  Date 

99 

Total 

538 

Total 

538 
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2.  There  is  no  systematic  way  of  analyzing  the  degree  of  uncer- 
tainty and  its  effect  on  the  safety  of  a design. 

3.  Additional  information  obtained  through  more  extensive  ex- 
plorations, better  quality  control  measures,  or  refined  loading  esti- 
mates cannot  be  used  in  evaluating  the  safety  of  a dam. 

In  a realistic  sense,  the  safety  of  a dam  is  a matter  of  acceptable 
risk,  and  within  this  premise,  the  logical  and  quantitative  basis  for 
evaluating  dam  safety  requires  the  concepts  and  methods  of  probability. 
The  safety  of  simple  engineering  structures  can  easily  be  determined  by 
comparing  the  resistance  of  the  structure  to  the  loads  or  load  effects 
applied  to  it.  Then  for  prescribed  or  assumed  distributions  of  the  load 
and  resistance  functions,  expressions  can  be  developed  relating  the 
safety  of  the  structure  to  the  properties  of  the  resistance  and  load 
distributions.  For  larger  engineering  structures  such  as  dams,  both  the 
resistance  and  the  loading  functions  will  depend  on  several  variables. 
For  example,  the  resistance  of  a dam  to  sliding  is  a function  of  the 
dead  load,  the  geometry  of  the  dam  profile,  uplift  forces,  compressive 
and  shearing  strengths  of  the  foundation  materials,  etc.  Moreover,  each 
of  these  variables  is  subject  to  an  estimation  error  or  uncertainty. 
Techniques  for  incorporating  these  uncertainties  in  the  design  equations 
must  be  formulated  before  it  will  be  possible  to  more  fully  evaluate  the 
safety  of  a dam  or  to  calculate  the  actual  or  mean  safety  factor  for  a 
specified  risk. 


Object i ve 

The  objectives  of  this  report  are:  (1)  to  develop  a probabilistic 
concept  for  gravity  dam  analysis  in  order  to  evaluate  dam  safety  in 
terms  of  the  various  sources  of  uncertainty  underlying  the  design  param- 
eters, and  (2)  to  assess  the  effects  of  various  uncertainties  on  dam 
safety. 


Scope 

The  probabilistic  dam  analysis  concept  described  in  this  document 
is  applicable  to  concrete  gravity  dams,  and  in  particular,  the  non- 
overflow sections  of  straight  concrete  dams  at  normal  operating  condi- 
tions or  maximum  operating  conditions  (i.e.,  normal  operating  condition 
plus  earthquake  effects).  The  loadings  considered  in  the  stability 
analysis  include  the  reservoir  and  tailwater,  the  weight  of  the  struc- 
ture, internal  hydrostatic  pressure  (uplift)  at  the  base  of  the  dam,  and 
earthquake  forces.  The  earthquake  forces  consist  of  the  inertial  forces 
caused  by  the  horizontal  acceleration  of  the  dam  itself  and  the  hydro- 
dynamic  forces  resulting  from  the  reaction  of  the  reservoir  water 
against  the  dam.  Secondary  loadings  have  been  neglected,  i.e.,  ice, 
wind,  wave,  and  earth  and  silt  pressures  have  not  been  included; 
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however,  if  these  loadings  are  deemed  to  be  significant,  the  technique 
can  be  modified  to  include  them.  Under  the  prescribed  loading  condi- 
tion, the  foundation  stresses  must  not  produce  a compressive  failure  of 
the  foundation  material,  and  the  sliding  resistance  of  the  dam  must  be 
adequate  to  prevent  a shear  failure.  For  the  case  at  hand,  both  of 
these  stability  requirements  are  evaluated  at  the  base  of  the  dam,  i.e., 
the  interface  oetween  the  concrete  dam  and  the  foundation  material. 
Moreover,  the  dam  is  not  embedded  or  anchored  to  the  foundation  mate- 
rial, and  the  plane  of  the  base  is  horizontal. 


2 CONCRETE  GRAVITY  DAMS 


General 

The  term  "concrete  gravity  dam"  refers  to  a solid  concrete  dam 
which  depends  primarily  on  its  own  weight  and  adhesion  with  the  founda- 
tion to  insure  stability  against  the  effects  of  all  imposed  forces. 
Gravity  dams  are  usually  straight,  but  may  be  slightly  curved  in  plan. 
The  general  shape  of  the  maximum  nonoverflow  section  of  a straight  grav- 
ity dam  is  shown  in  Figure  1.  The  shape  of  the  section  is  determined  by 
the  prescribed  loading  conditions,  the  shear  resistance  of  the  rock,  and 
the  height  of  the  section.  The  upstream  face  of  the  gravity  dam  is  usu- 
ally vertical  to  concentrate  the  concrete  in  a location  where  it  acts  to 
overcome  the  effects  of  the  reservoir  water  load.  The  lower  portion  of 
the  upstream  face  may  be  battered  to  increase  the  thickness  at  the  base 
in  order  to  improve  the  base's  sliding  safety.  The  downstream  face  will 
usually  have  a uniform  slope  which  is  determined  by  both  the  stress  and 
stability  requirements  at  the  base.  The  crest  thickness  may  be  dictated 
by  roadway  or  other  requirements,  such  as  possible  ice  pressures  or 
impact  from  floating  objects.  Generally,  the  downstream  face  is  verti- 
cal from  the  downstream  edge  of  the  crest  to  the  intersection  with  the 
sloping  downstream  face. 
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Design  Forces 


Several  forces  must  be  considered  in  the  design  of  gravity  dams. 

The  nature  of  most  of  these  forces,  unfortunately,  is  such  that  an  exact 
determination  cannot  be  made.  Thus,  the  magnitude,  direction,  and 
location  of  these  forces  must  be  assumed  after  consideration  of  all  the 
available  facts  and,  to  a certain  extent,  must  be  based  on  judgment  and 
experience.  The  principal  forces  commonly  considered  in  the  design  of  a 
concrete  gravity  dam  are  the  following. 


1.  Dead  Load.  The  dead  load  is  the  weight  of  concrete  plus  appur- 
tenances such  as  gates  and  bridges.  In  the  determination  of  the  dead 
load,  relatively  small  voids,  such  as  galleries,  normally  are  not 
deducted.  The  unit  weight  of  concrete  is  generally  assumed  to  be  150 
pcf  unless  specific  data  are  available  from  a concrete  mix  analysis. 

2.  Reservoir  and  Tailwater  Pressures.  Reservoir  and  tailwater 
loads  are  obtained  from  reservoir  operations  studies.  A triangular  dis- 
tribution of  the  static  water  pressure  acting  normal  to  the  face  of  the 
dam  is  assumed  in  the  design  of  most  dams.  Although  the  weight  of  water 
varies  slightly  with  temperature,  the  weight  usually  adopted  in  the 
design  of  dams  is  62.5  pcf. 
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Figure  1.  General  profile  of  nonoverflow  section  of  a gravity  dam. 

3.  Internal  Hydrostatic  Pressure.  Dams  are  subject  to  water  pres- 
sure not  only  on  exposed  faces  but  also  on  their  bases.  These  internal 
pressures  produce  uplift,  which  causes  a reduction  in  the  effective 
weight  of  the  structure.  The  distribution  of  the  internal  hydrostatic 
pressure  along  a horizontal  section  through  a gravity  dam  is  usually  as- 
sumed to  vary  linearly  from  full  reservoir  pressure  at  the  upstream  face 
to  zero  or  tailwater  pressure  of  the  downstream  face  and  to  act  over  the 
entire  area  of  the  section.  The  pressure  distribution  may  be  adjusted 
to  reflect  the  size,  location,  and  spacing  of  drains,  if  they  are  used. 
When  the  line  of  drains  intersects  the  foundation  near  the  upstream 
face,  which  is  consistent  with  the  requirements  for  gallery  locations, 
the  pressure  distribution  is  normal ly  assumed  to  vary  linearly,  but  is 
adjusted  by  a hydrostatic  pressure  intensity  factor,  K,  as  shown  in 
Figure  2. 

4.  Earthquake  Forces.  In  regions  where  earthquakes  occur,  dams 
must  resist  the  inertia  effects  caused  by  the  sudden  movement  of  the 
earth’s  crust.  Two  types  of  earthquake  forces  are  produced:  the 

^ EM  1110-2-2200,  Gravity  Dam  Design  (U.S.  Army  Corps  of  Engineers, 

1%8). 


! 


Figure  2.  Distribution  of  internal  hydrostatic  pressure.  From  EM  1110-2-2200, 
Gravity  Dam  Design  (U.S.  Army  Corps  of  Engineers,  1958). 


inertial  force  due  to  the  acceleration  of  the  mass  of  the  dam  and  the 
hydrodynamic  force  resulting  from  the  reaction  of  the  reservoir  water  on 
the  dam.  For  a gravity  dam  with  a full  reservoir,  the  most  unfavorable 
direction  of  the  earthquake  is  upstream  normal  to  the  axis,  which  pro- 
duces horizontal  forces  which  act  downstream. 

The  force  imposed  on  the  mass  of  the  dam  may  be  determined  by  the 
equation 

PQ  = a WQ  [Eq  1]  1 

where  PQ  = inertial  force  due  to  earthquake  per  unit  length 
a * ratio  of  the  earthquake  acceleration  to  gravity 

; 

W0  = total  dead  weight  of  the  dam  per  unit  length 


The  increase  in  water  pressure  on  the  dam  caused  by  the  earthquake 
acts  concurrently  with  the  above  force.  The  distribution  of  these 
increased  pressures  is  usually  assumed  to  be  parabolic,  and  the  force 
may  be  computed  using  the  following  equation 
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Pe  = 2/3  Co  [Eq  2] 

where  Pp  = total  hydrodynamic  force  due  to  earthquake  per  unit  length 
= height  of  the  reservoir  water 

C = factor  depending  principally  on  the  height  of  the 

dam  and  the  earthquake  period,  t , as  shown  in  Figure  3. 

Typically,  an  acceleration  of  O.lg,  i.e.,  u = 0.1,  has  been  used  in 
the  regions  where  moderate  to  severe  shocks  have  been  experienced  or 
where  conditions  are  unfavorable.  In  more  favorable  regions  this  factor 
has  been  reduced  to  0.05g.  The  period  of  vibration  of  the  earthquake, 
t , is  usually  assumed  to  be  1 second. 

5.  Foundation  Reactions.  In  general,  the  resultant  of  all  hori- 
zontal and  vertical  forces  on  the  dam  including  uplift,  must  be  balanced 
by  an  equal  and  opposite  reaction  at  the  foundation,  consisting  of  the 
total  vertical  reaction  and  the  total  horizontal  shear  and  friction. 
Moreover,  for  the  dam  to  be  in  static  equilibrium,  the  location  of  this 
force  must  be  such  that  the  summation  of  moments  is  equal  to  zero.  The 
distribution  of  the  vertical  reaction  is  assumed  to  be  trapezoidal  for 
convenience  only,  with  the  knowledge  that  the  elastic  and  plastic  prop- 
erties of  the  foundation  material  and  the  concrete  affect  the  actual 
distribution.  The  problem  of  determining  the  actual  distribution  is 
complicated  by  the  horizontal  reaction,  internal  stress  relations,  and 
other  theoretical  considerations.  Also,  the  variation  of  the  foundation 
materials  with  depth,  and  geological  discontinuities,  which  disrupt  the 
tensile  and  shearing  resistance  of  the  foundation,  make  the  problem  more 
complex.  The  maximum  foundation  pressure  at  any  point  is  equal  to  the 
computed  foundation  reaction  plus  the  uplift  pressure. 

In  addition  to  the  principal  forces  described  above,  there  are  sec- 
ondary forces  that  may  be  considered  where  warranted.  These  include  the 
fol  1 owi ng : 

1.  Earth  and  Silt  Pressures.  In  the  design  of  abutment  sections, 
some  depths  of  backfill  usually  must  be  considered.  Silt  pressures  on 
both  the  upstream  and  downstream  faces  of  the  dam  may  develop  during  the 
expected  useful  life  of  the  dam  and  should  be  given  some  consideration. 

2.  Ice  Pressure.  If  ice  will  develop  to  appreciable  thickness  and 
remain  for  a long  duration,  ice  pressure  should  be  considered  in  the 
design  of  the  dam.  For  dams  in  this  country,  ice  pressures  usually  will 
not  exceed  10,000  plf  of  the  dam  and  should  be  applied  at  normal  pool 
elevation. 
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Valuts  of  C 


Figure  3.  Value  of  earthquake  factor  C. 


3.  Wave  Pressures.  In  some  instances,  wave  pressures  can  have  an 
effect  on  the  dam.  However,  usually  the  height  of  the  waves  is  a more 
important  factor  since  it  establishes  the  freeboard  requirement  for  the 
dam. 


4.  Wind  Pressures.  Wind  pressures  of  30  psf  are  usually  used  in 
the  design  of  dams.  However,  unless  there  are  large  appurtenances  or 
freeboard  requirements,  these  forces  are  small. 


Failure  Mechanisms 

There  are  two  direct  ways  in  which  a gravity  dam  may  fail: 

1.  By  sliding  on  a horizontal  or  nearly  horizontal  joint  above  the 
foundation,  on  the  foundation,  or  on  a horizontal  or  nearly  horizontal 
seam  in  the  foundation. 

2.  By  overturning  on  a horizontal  joint  within  the  dam,  at  the 
base,  or  at  a plane  below  the  base. 
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The  direct  cause  of  sliding  is  the  presence  of  horizontal  forces 
greater  than  the  combined  shearing  resistance  of  the  joint  or  base  and 
the  static  friction  induced  by  the  vertical  forces. 

The  direct  cause  of  pure  overturning,  not  preceded  by  some  other 
type  of  failure,  if  tension  is  ignored,  is  the  presence  of  horizontal 
forces  great  enough  compared  to  the  vertical  forces  to  cause  the  re- 
sultant of  all  forces  acting  on  the  dam  above  any  horizontal  plane,  in- 
cluding uplift,  to  lie  outside  the  limits  of  the  dam. 

As  the  resultant  approaches  the  face,  the  compressive  stresses 
increase  rapidly,  hence  overturning  would  be  preceded  and  accelerated  by 
a compression  failure.  In  fact,  a dam  with  the  resultant  well  inside 
the  base  may  overturn  if  the  toe  of  the  dam  fails  due  to  crushing  or 
other  causes  and  reduces  the  effective  length  of  the  joint  or  base  suf- 
ficiently to  cause  the  resultant  to  lie  outside  the  effective  base. 

A dam  may  start  to  overturn  but  finally  fail  by  sliding.  If  the 
resultant  passes  appreciably  outside  the  middle  third  (or  kern,  if  the 
base  is  of  irregular  form),  a horizontal  tension  crack  may  occur,  which 
reduces  the  shearing  strength  of  the  joint  or  base.  Also,  the  admission 
of  reservoir  water  pressure  to  the  fissure  increases  the  uplift,  reduc- 
ing the  net  reaction  and  the  frictional  resistance  to  horizontal  motion. 
As  a result,  sliding  may  occur. 


Loading  Conditions 

The  following  loading  conditions  and  requirements  are  generally 
considered  in  the  design  of  moderate  height  gravity  dams;3 

1.  Construction  Condition.  Dam  completed  but  no  water  in  reser- 
voir, no  tailwater,  wind  load  on  downstream  face. 

2.  Normal  Operating  Condition.  Pool  elevation  at  top  of  closed 
spillway  gates,  where  spillway  is  gated,  and  at  spillway  crest,  where 
spillway  is  ungated.  Minimum  tailwater  elevation  for  gated  and  ungated 
spillways.  Ice  pressure,  if  applicable. 

3.  Induced  Surcharge  Condition.  Pool  elevation  at  top  of  par- 
tially opened  gate.  Tailwater  pressure  at  full  value  for  nonoverflow 
section  and  60  percent  of  full  value  for  overflow  section.  Ice  pres- 
sure, if  appl icable. 


1110-2-2200,  Gravity  Dam  Design  (U.S.  Army  Corps  of  Engineers, 
1958). 


16 


1 


item 


■m 


z:\ 


4.  Flood  Discharge  Condition.  Reservoir  at  maximum  flood  pool 
elevation.  All  gates  open  and  tailwater  at  flood  elevation.  Tailwater 
pressure  at  full  value  for  nonoverflow  section;  60  percent  of  full  value 
for  overflow  sections  for  all  conditions  of  deep  flow  over  spillway, 
except  that  full  value  will  be  used  in  all  cases  to  compute  the  uplift. 
No  ice  pressure. 

5.  Construction  Condition  With  Earthquake.  Earthquake  acceler- 
ation in  a downstream  direction.  No  water  in  reservoir.  No  wind  load. 
No  tailwater. 

6.  Normal  Operating  Condition  With  Earthquake.  Earthquake  accel- 
eration in  an  upstream  direction.  Reservoir  at  top  of  closed  gate  for 
gated  spillways  and  at  spillway  crest  for  ungated  spillways.  Minimum 
tailwater.  No  ice  pressure. 


Stability  Requirements 

The  basic  requirements  for  stability  of  a gravity  dam  for  all  con- 
ditions of  loading  are: 

1.  That  it  be  safe  against  overturning  at  any  horizontal  plane 
within  the  dam,  at  the  base,  or  at  a plane  below  the  base. 

2.  That  it  be  safe  against  sliding  on  any  horizontal  plane  within 
the  dam,  on  the  foundation,  or  on  any  horizontal  or  nearly  horizontal 
seam  in  the  foundation. 

3.  That  allowable  unit  stresses  in  the  concrete  or  in  the  founda- 
tion material  shall  not  be  exceeded. 


Sliding  and  Overturning  Criteria 

The  horizontal  forces  acting  on  a dam  tend  to  cause  the  structure 
to  slide  along  horizontal  planes.  Resistance  to  sliding  at  or  below  the 
base  is  a function  of  the  shearing  strength  of  the  foundation,  or  at  any 
plane  above  the  base,  the  shearing  strength  of  the  concrete  in  the  dam. 
Experience  has  shown  that  the  shearing  resistance  of  the  foundation  or 
concrete  need  not  be  investigated  if  the  ratio  of  horizontal  forces  to 
vertical  forces  (ZH/Z/V)  is  such  that  a reasonable  safety  factor  re- 
sults. This  will  require  that  the  ratio  of  EH/EV  be  well  below  the  co- 
efficient of  sliding  friction  of  the  material. 

When  the  resultant  of  all  forces  acting  above  any  horizontal  plane 
through  a dam  intersects  that  plane  outside  of  its  middle  third  or  kern, 
tensile  stresses  are  set  up.  Since  the  tensile  strength  of  concrete  is 
usually  considered  to  be  zero,  it  is  generally  required  for  all  normal 
loading  that  the  resultant  intersect  the  plane  within  the  middle  third 


or  kern.  However,  because  of  the  short  duration  of  earthquake  forces, 
it  is  considered  justifiable  to  modify  this  requirement  when  these 
forces  are  included. 

Table  2 indicates  acceptable  limits  of  the  location  of  the  re- 
sultant and  the  maximum  sliding  factor  for  all  loading  conditions. 
Where  the  sliding  factor  exceeds  the  values  given,  the  shearing  re- 
sistance should  be  investigated. 


Shear  Friction  Factor  of  Safety 

Where  the  ratio  of  T.H/>;V  indicates  the  need  to  investigate  the 
shearing  resistance  to  horizontal  movement,  the  shear-friction  formula 
should  be  used  to  determine  the  factor  of  safety  inherent  in  the  founda- 
tion or  concrete.  This  factor  of  safety  may  be  evaluated  as  follows: 


/ <-  i _ _Q.  IjV  _taru{>  + sA 

^ s-f;  " TH  ' i:h 


[Eq  3] 


where  (S  = 
s = 
A = 

Q = 
= 


shear-friction  safety  factor 
unit  shearing  strength  of  material 
area  of  the  horizontal  plane  considered 
sliding  resistance  of  the  dam  per  unit  length 
angle  of  internal  friction  of  foundation  material 


A minimum  shear-friction  safety  factor  (S  of  4 is  generally  required 
for  the  maximum  loading  condition. 


Analysis  Method 

Generally,  the  Gravity  Method  of  stress  and  stability  analysis  is 
used  in  the  preliminary  studies  of  gravity  dams  and  the  final  designs  of 
straight  gravity  dams  in  which  the  transverse  contraction  joints  are 
neither  keyed  nor  grouted.  The  gravity  method  provides  an  approximate 
determination  of  stresses  in  a section  of  a gravity  dam.  Each  section 
of  the  dam  is  assumed  to  act  independently  and  is  treated  as  if  it  were 
a cantilever  with  the  loads  transmitted  to  the  foundation  by  cantilever 
action  and  resisted  by  the  weight  of  the  cantilever.  The  forces  acting 
on  the  cantilever  element,  including  uplift  and  earthquake,  are  shown 
for  the  normal  loading  conditions  in  Figure  4.  Symbols  used  to  define 
this  loading  condition  are  presented  in  the  appendix. 
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Limits  for  the  Location  of  the  Resultant  and  Maximum  Sliding  Factor 
(From  LM  1110-2200,  Gravity  Dam  Design  (U.S.  Army  Corps  of  Engineers,  1958) 

Maximum 


..oading  Condition 

Resultant  to  Fall  Within  Following 
Central  Percent  of  Base  Width 

SI iding 
Factor  EH/EV 

I 

33  1/3 

_ „ _ _ 

II 

33  1/3 

0.65 

III 

33  1/3 

0.65 

IV 

33  1/3 

0.65 

V 

(For  conditions  V and  VI,  within 

— 

VI 

base  but  allowable  foundation 
pressure  must  not  be  exceeded) 

0.85 

Forces  and  Moments 

The  resultants  of  the  vertical  and  horizontal  forces  acting  on  the 
section  of  the  dam  shown  in  Figure  4 can  be  readily  computed.  Assuming 
positive  forces  and  moments  are  as  indicated  by  the  directional  arrows 
shown  in  Figure  4,  the  equations  for  the  resultant  of  the  vertical  and 
horizontal  forces  are,  respectively. 


M I ■ 


Mq=  4/15CaH^  + 1/6^  - 1/6ywh|  + otYc ( 1 / 6L, h|+1  / 2L2H^+1  /6L3H^ ) 


-1/2>wL1(H1-H3)(L-L1)  - 1/2ywL1H3(1/2L-1/3L1  ) - 1 / 2y CL H3 (1  / 2L- 2/ 3L 3 ) 


-YCL2H0(1/2L-Lr1/2L2)  + 1/2ycL3H40/2L-2/3L3)  + l/Zy^tane' 


(l/2L-l/3H2tane’ ) + 1 / 1 2ywK(H] -H2 )L" 


[Eq  6] 


In  turn,  the  vertical  foundation  pressure  can  be  determined  at  any 
point  on  the  base  from  the  equations 


p,=ZV  0 


v 

I 


[Eq  7] 


p - ?V  + V 
r A 


[Eq  8] 


where  PJ.,  P"  = unit  vertical  reaction  in  the  foundations  at  the 

upstream  and  downstream  faces  of  the  dam,  respectively 


m = distance  from  the  center  of  gravity  of  the  base  to  the 
point  under  consideration. 


For  a section  of  a gravity  dam  of  unit  length,  certain  sim- 
plifications are  possible.  The  base  of  the  section  is  a rectangle  of 


unit  width  and  length  L;  therefore,  A = L.  The  center  of  gravity  is  at 
, and  m = L/2,  and  I = 1/12  l3.  Substituting  these  values 


the  midpoint 
in  the  above  expression  yields 


, _ £V 


P'  = 
r L 


6M 

o 

T*" 


[Eq  9] 


p - SV  + 6M0 

Pr  L "^T 


[Eq  10] 
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To  obtain  the  maximum  vertical  foundation  pressures,  the  uplift 
pressures  must  be  added  to  the  vertical  foundation  pressure. 


P' 

v 


r\ 

L 


+ p1 
u 


[Eq  11] 


P" 

v 


P" 

u 


[Eq  12] 


where  P^,  P^  = maximum  vertical  unit  pressure  on  the  upstream  and  down- 
stream faces  of  the  dam,  respectively 

Pu’  Pu  = unit  uPlift  pressures  at  the  base  of  the  dam  on  the  up- 
stream and  downstream  faces  of  the  dam,  respectively. 


The  maximum  vertical  compressive  stresses,  however,  are  not  the 
maximum  stresses  which  occur  in  the  dam.  The  maximum  stresses  occur  at 
the  ends  of  joints,  on  inclined  planes,  normal  to  the  face  of  the  dam. 

In  the  general  case  where  external  water  forces  are  involved  due  to  the 
reservoir  and  tailwater,  the  corresponding  inclined  pressures  are 

PI  = P^secV  - P^tan20'  [Eq  13] 

and 


P-  = PySecV  - P^tane"  [Eq  14] 


where  PI,  PV  = included  pressure  at  the  upstream  and  downstream  faces, 
respectively 

O',  O"  = angle  between  the  upstream  and  downstream  faces  and  the 
vertical,  respectively 

Pn’  Pri  = external  normal  pressures  at  the  upstream  and  downstream 
faces,  respectively. 
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Whether  the  maximum  pressure  in  the  foundation  is  equal  to  the 
inclined  pressure  in  the  dam  at  the  base  or  to  the  maximum  vertical 
foundation  pressure  is  frequently  the  subject  of  debate.  Undoubtedly, 
stress  conditions  change  rapidly  in  the  rock  beneath  the  toe  of  the  dam, 
but  it  seems  rational  and  on  the  side  of  safety  to  assume  that  the 
maximum  stress  in  the  rock  in  immediate  contact  with  the  base  of  the  dam 
equals  the  inclined  toe  stress  in  the  dam. 


3 RELIABILITY  THEORY 


General 

In  the  classical  concept  of  reliability,  the  random  resistance,  R, 
and  the  random  load  (or  load  effect),  S,  are  known  and  the  uncertainties 
associated  with  the  inherent  randomness  are  described  by  known  proba- 
bility distributions.  Eor  the  practical  case  at  hand,  this  concept  is 
inadequate.  First,  R and  S are  not  known,  and  it  is  only  possible  to 
predict  these  variables  using  theoretical  models  R and  S,  respectively. 
Such  predictions  are  invariably  imperfect,  thus  introducing  additional 
uncertainties.  Second,  the  uncertainties  underlying  the  resistance  and 
load  effects  must  be  given  in  terms  of  quantitative  measures. 


Extended  Reliability  Theory 

Recent  developments  in  reliability  theory  have  extended  the  classi- 
cal concept  to  yield  a first-order  approximation  for  practical  applica- 
tions, such  as  computing  the  probability  of  failure  of  a concrete  grav- 
ity dam.  4 The  essence  of  this  theory  and  a synopsis  of  the  relevant 
equations  for  computing  the  means,  uncertainties,  and  probability  of 
failure  are  presented  below. 

1.  Mean  and  Uncertainty.  Let  Y be  a random  variable  representing, 
for  example,  the  resistance,  R,  or  the  applied  load  effect,  S,  in  a 
structure.  Invariably,  Y is  a function  of  other  variables;  e.g., 

Y = f(Xr  X2,...,  Xn)  [Eq  15] 


Presumably,  the  model  for  Y,  as  represented  by  the  function  in  Eq 
15,  as  well  as  Xi,...,Xn,  would  represent  reality  exactly.  In 
engineering  practice,  however,  this  generally  is  not  the  case;  f and 
Xl,...,Xn  must  be  predicted  or  estimated,  and  thus  are  subject  to  pre- 
diction errors.  To  adjust  for  an  imperfect  prediction,  corrective  fac- 
tors Nf  and  are  introduced,  such  that 


f - 


Nf  f 


[Eq  16] 


4 A.  H-S.  Ang  and  C.  A.  Cornell,  "Reliability  Bases  of  Structural  Safety  and 
Desiqn,"  ASCE,  Journal  of  the  Structural  Division,  Vol  100,  No.  St  9 
(1974),  pp~T755^l 769;  
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Xi  = NX.  Xi 


[Eq  17] 


where  fjs  the  empirical  or  theoretical  function  adopted  as  a model  of 
f,  and  Xj  is  the  mode^  of  Xj . _Here,  Nf,  NXj , and  Xj  are  random  vari- 
ables with  means  Nf,  • 'X-j , and  X-j , and  coefficients  of  variation  (c.o.v.) 
<5f>  AXj,  and  “X-j,  respectively.  The  uncertainties  associated  with  the 
basic  variability  in  X-j  are,  therefore,  measures  by  °Xj , whereas  AXj 
represents  the  prediction  uncertainty  in  Xj . Consistent  with  the  first- 
order  approximation,  AXj  will  be  ascribed  entirely  to  the  uncertainty  in 
the  predicted  Xj.  Furthermore,  the  mean  values  Nf  and  ftyj  represent,  re- 
spectively, the  bias  in_the  model  f (usually  a deterministic  function) 
and  the  estimated  mean  X . . 


By  first-order  approximation,  the  total  c.o.v.  of  X,  then  (by 
virtue  of  Eq  17)  is 


\ =A2. + Ax2i 


[Eq  18] 


Similarly,  the  total  c.o.v.  of  f is 


= A2  . .2 
f / 6^  + 


[Eq  19] 


where  5f  represents  the  basic  variability  about  the  model  function  f;  and 
Af  denotes  any  imperfections  in  the  form  of  the  equation  used. 

Substituting  Eqs  16  and  17  into  Eq  15,  and  using  a first-order  ap- 
proximation, the  mean  and  c.o.v.  of  Y are  easily  found  to  be 


Y - Nf  ftN^*,,  Y2 


„2  ♦ 5 r?  A2  if2  x2  a2 

f y2  i£1  9X..  NX.  Ai  “X. 


[Eq  20] 


[Eq  21] 


+ 1 l ( ax'.)  (axf .)  px  ’X.  Nx.Nx.xixjux.^x. 

j=l  j=l  1 J 1 J 1 J 1 J 


. ...  . ..  . 'HwU 


where  the  derivations  are  evaluated  at  the  mean  value  and  . u is  the 


correlation  coefficient  between  Xj  and  X^ 


V «j 


2.  Probability  of  Failure.  The  probability  of  failure  is  defined  as 

r\ 


n 


2 


-tn 


/ 2 2 
V WR  + S2s 


[Eq  22] 


where  F = cumulative  distribution  function  of  n9 
n2  d 


!n(s)~  4s) 


Ur  + U$ 


R,S  r true  random  resistance  and  true  random  load,  respectively 

a /s 

R,S  = models  of  R and  S 

R,S  = mean  values  of  R and  S,  respectively 


Ur  = /6R+AR  = coefficient  of  variation,  representing  measure  of  the  total 
uncertainty  of  resistance,  R 

°R 

6R  = — - = coefficient  of  variation,  rgpresenting  measure  of  the  basic 
R variability  of  resistance,  R 

AR  = coefficient  of  variation  of  resistance,  representing  a 
measure  of  the  errors  in  predicting  I? 

/2  2 

s*S  r •'°s+As  - coefficient  of  variation,  representing  measure  of  the  total 
uncertainty  of  load,  S 


6c  = ~ = coefficient  of  variation,  representing  measure  of  the  basic 
S variability  of  load,  S 

Ac  * coefficient  of  variation  of  load,  representing  a measure  of 
the  errors  in  predicting  S. 
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Eq  22  is  to  be  used  when  the  total  uncertainties  are  not  large  (j<0.3). 
Otherwise,  a more  general  form  is 


/'  ♦«! 


/,T 


(l  + fip)(i  + n*) 


CEq  23] 


Distribution  Function 


The  probability  of  failure,  as  given  in  Eq  22,  depends  on  the  dis- 
tribution functions  for  R and  S.  Consequently,  the  probability  of  fail- 
ure may  be  sensitive  to  the  assumed  distribution  functions  if  the  risk 
is  very  small  (e.g.,  < 10'5).  This  is  a consequence  of  the  behavior  of 
the  extreme  tails  of  the  distribution  functions.  Unfortunately,  when 
the  available  data  are  limited  to  the  central  range  of  the  variates,  it 
is  virtually  impossible  to  determine  the  validity  of  one  distribution 
function  versus  several  others.  Moreover,  the  choice  of  the  distribu- 
tion should  recognize  the  practical  considerations  such  as  convenience 
of  implementation  and  availability  of  probabilistic  information.  In 
similar  engineering  appl ications,  the  log-normal  distribution  has  been 
favored;  therefore,  it  was  selected  for  this  study. 
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L\  PROBABILISTIC  LOAD  AND  RESISTANCE  MODELS 


General 

In  the  stability  analysis  of  concrete  gravity  dams,  the  horizontal 
force,  sliding  resistance,  and  foundation  pressure  are  functions  of 
several  variables.  For  example,  the  horizontal  force  is  a function  of 
the  reservoir  and  tailwater  elevations,  the  earthquake  acceleration  co- 
efficient, and  the  earthquake  factor.  The  sliding  resistance  is  a func- 
tion of  the  reservoir  and  tailwater  elevations,  the  internal  hydrostatic 
pressure  (i.e.,  uplift),  the  angle  of  internal  friction  of  the  founda- 
tion material,  and  the  unit  shear  stress  (cohesion)  along  the  failure 
plane.  Likewise,  the  foundation  pressure  is  a function  of  many  of  these 

variables.  Since  the  profile  of  a nonoverflow  section  of  a dam  is  usu- 

ally well  defined,  the  dimensions  of  the  section  will  be  considered  con- 
stant. Likewise,  the  unit  weights  of  concrete  and  water  will  be  consid- 
ered constant,  since  their  basic  variabilities  are  quite  small. 
Therefore,  the  principal  random  variables  in  the  stability  analysis  of 
concrete  gravity  dams  are: 

C = the  earthquake  factor 

= height  of  the  reservoir  water 

= height  of  the  tailwater 

K = hydrostatic  pressure  intensity  factor 

s = unit  shear  strength  along  the  failure  plane 

o = unit  compressive  strength  of  the  foundation  material 

« = earthquake  acceleration  coefficient 

<f>  = angle  of  internal  friction  of  the  foundation  material 

To  apply  the  extended  reliability  theory  formulated  in  the  previous 
section  to  concrete  gravity  dams,  probabilistic  models  of  the  load  and 
resistance  functions  for  the  sliding  and  overturning  modes  of  failure 
must  be  developed  which  incorporate  corrective  factors  and  coefficients 
of  variation  terms  to  account  for  basic  variabilities  and  prediction  un- 
certainties. Moreover,  some  of  the  sources  for  the  basic  variabilities 
and  prediction  uncertainties  associated  with  the  principal  variables 
need  to  be  identified,  although  generally  they  cannot  be  quantified 
until  specific  data  and  information  for  an  actual  dam  are  available. 
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Horizontal  Force  Model 


Equations  for  Mean  and  Uncertainty 

The  horizontal  forces  tending  to  cause  sliding  of  the  dam  are  sum- 
marized in  Eq  5.  The  major  variables  in  this  equation  are  H],  C,  u,  and 
H2 ; thus,  correction  factors  , Nq,  N(1,  nH?,  respectively,  and  must 
be  introduced  into  Eq  5 to  yield  the  mean  value  of  total  horizontal 
force,  H,  as 


H - Nh[2/3NaNcN^aCH2  * [Eq  24] 

Here  and  in  succeeding  expressions,  bars  above  the  random  variables  to 
denote  mean  values  have  been  deleted  to  simplify  the  presentation;  how- 
ever, it  should  be  understood  that  mean  values  of  these  random  variables 
are  implied. 

To  evaluate  the  uncertainty  associated  with  the  total  horizontal 
force,  it  is  necessary  to  evaluate  Eq  21.  For  the  purpose  of  this 
study,  H-j  and  H2  are  assumed  to  be  correlated  and  the  correlation  coef- 
ficients between  the  other  variables  are  assumed  to  be  zero;  i.e.,  they 
are  statistically  independent.  Thus,  the  total  uncertainty  associated 
with  the  total  horizontal  force  is 

N2 

«H  ■ ah  * 7»/3“CHl+\"l>2  "h^H,  + hwHZ>2  NH2H2!1H2 

+ (2/3CH2+Wq)2  N2aV  + 2/3(aH2)2  N2C2^2  [Eq  25] 


- (4/3aCHi+YwH1)(YwH2)pHiH^NH^NH^H1H2ilH^H_  ] 


Bias  and  Uncertainty  of  Principal  Variables 

Height  of  Reservoir  and  Tailwater.  Most  gravity  dams  will  have 
gated  spillways  or  outlet  works  to  control  the  discharge  of  the  reser- 
voir water  for  irrigation,  electrical  power  generation,  downstream  water 
control,  etc.  Normally,  there  will  not  be  large  fluctuations  in  the 
reservoir  and  tailwater  elevations.  Consequently,  the  bias  and  uncer- 
tainty associated  with  the  reservoir  and  tailwater  elevations  will  be 
quite  small,  though  finite. 
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Earthquake  Acceleration  Coefficient.  In  the  pseudo-dynamic  ap- 
pro a cTTTo  eaHTiquaireTbli^^  gravity  dams,  additional  static 

forces  are  included  in  the  basic  design  equations  to  represent  the 
inertial  forces  on  the  body  of  the  dam  and  the  water  in  the  reservoir 
during  an  earthquake.  Typically,  an  earthquake  acceleration  coefficient 
of  0.05  g to  0.1  g has  been  adopted  to  account  for  the  effects  of  hori- 
zontal accelerations,  and  the  effects  of  vertical  accelerations  are  fre- 
quently ignored.  Recently,  tne  seismic  zoning  maps  for  the  United 
States  were  updated,  and  the  earthquake  coefficients  in  a number  of 
areas  were  increased  to  a maximum  value  of  0.15  g (Figures  5 through 
7.) ^ However,  the  problem  remains  that  these  maps  are  based  on  estimates 
of  the  maximum  ground  shaking  experienced  during  the  recorded  historical 
period  with  only  minor  consideration  of  how  frequently  such  motions  have 
occurred . 

What  is  actually  needed  in  planning  the  design  of  a dam  in  a region 
of  potential  earthquake  activity  is  a specification  of  the  maximum 
ground  motion  expected  at  the  site  during  the  useful  life  of  the  dam; 
i.e.,  a site-dependent  seismic  hazard  analysis  should  be  conducted  to 
synthesize  all  regional  information  and  general  knowledge  about 
earthquakes  to  obtain  the  probabilities  associated  with  the  occurrence 
of  different  motions  at  the  site.  Such  analysis  should  realistically 
model  the  fault  system,  the  occurrence  rate  for  earthquake  magnitudes  of 
interest,  and  the  attenuation  relationship  from  the  source  to  the  site; 
in  addition,  the  analysis  should  account  for  uncertainty  due  to  the 
degree  of  scatter  in  the  data  and  the  prescribed  relationships  used  in 
the  model.  Through  such  analysis,  the  annual  probabilities  (or  return 
periods)  associated  with  various  intensities  of  ground  motion  at  the 
site  can  be  derived  systematically  and  objectively,  using  the  seismic 
information  available  for  a region  and  the  relevant  physical  re- 
lationships of  earthquakes.  Moreover,  by  computing  the  hazard  at  vari- 
ous sites  in  the  region,  seismic  hazard  maps  can  be  constructed  in  the 
form  of  isoseismal  contours  corresponding  to  given  time  periods  and 
probabilities  of  non-exceedance. 

Recently,  the  Applied  Technology  Council,  an  affiliate  of  the 
Structural  Engineers  Association  of  California,  made  a first  attempt  at 
incorporating  seismic  risk  procedures  into  zoning  maps  for  the  United 
States  as  part  of  the  Cooperative  Federal  Program  in  Building  Practices 
for  Disaster  Mitigation.  These  seismic  hazard  maps  were  developed  to 
incorporate  the  attenuation  of  ground  motions  with  distance  to  the  site 
from  the  anticipated  earthquake  sources  and  to  have  the  probability  of 
exceeding  the  design  ground  shaking  roughly  the  same  in  all  parts  of  the 
country.  The  maps  did  not,  however,  include  detailed  representation  of 
the  fault  system  or  provisions  for  incorporating  uncertainty.  To  re- 
present the  intensity  of  ground  shaking,  two  parameters  were  adopted: 


ER  1110-2-1806,  Earthquake  Design  and  Analysis  for  Corps  of  Engineer 
Dams  (U.S.  Army  Corps  of  Engineers,  1977). 
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Figure  5.  Seismic  zone  map  for  contiguous  states.  From  ER  1110-2-1806,  Earthquake  Design  and 
Analysis  for  Corps  of  Engineers  Dams  (U.S.  Army  Corps  of  Engineers,  1977). 
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Figure  7.  Seismic  zone  map  for  Hawaii.  From  ER  1110-2-1806,  Earthquake  Design 
and  Analysis  for  Corps  of  Engineers  Dams  (U.S.  Army  Corps  of  Engineers 


effective  peak  acceleration  (EPA)  and  effective  peak  velocity-related 
acceleration  (EPV).  The  maps  of  EPA  and  EPV  are  shown  in  Figures  8 and 
9,  respectively. 

In  Figure  8,  the  maximum  value  assigned  to  EPA  anywhere  on  the  map 
is  0.4  y;  however,  this  maximum  is  not  intended  to  represent  an  absolute 
upper  bound  of  EPA.  There  are  locations  inside  the  0.4  g contour  where 
higher  values  of  EPA  would  be  appropriate;  however,  contouring  such 
small  areas  would  have  required  special  site-dependent  studies  similar 
to  those  previously  mentioned. 

The  map  of  EPV  shown  in  Figure  9 was  constructed  by  modifying  the 
EPA  map;  however,  a different  attenuation  relationship  was  used  in  the 
eastern  half  of  the  country.  In  the  western  half  of  the  country,  the 
distance  required  for  the  EPV  to  reduce  to  one-half  the  original  value 
from  a large  earthquake  is  about  80  miles.  Data  on  attenuation  of  Mod- 
ified Mercalli  Intensity  in  the  eastern  United  States  have  indicated 
that  within  100  miles  of  a large  earthquake,  attenuation  was  the  same  as 
in  the  west.  Thereafter,  the  distance  required  for  the  velocity  to  be 
reduced  by  one  half  almost  doubles.  These  modifications  resulted  in  the 
map  shown  in  Figure  9.  In  using  acceleration  and  velocity  coefficients 
from  these  maps,  it  must  be  remembered  that  forces  computed  using  these 
coefficients  are  conditional  upon  the  occurrence  of  an  earthquake  of 
sufficient  magnitude  to  produce  these  design  accelerations  or  velocities 
at  the  site.  This  conditional  probability  is  about  0.1  to  0.2. 

In  view  of  the  previous  discussion  and  the  disparity  between  Figure 
5 and  Figures  8 and  9,  the  earthquake  acceleration  coefficient  is  sub- 
ject to  considerable  bias  and  uncertainty. 

Earthquake  Factor.  In  recent  years,  a great  deal  of  research  has 
focused  on  the  interaction  between  the  dam  and  the  reservoir  water 
during  earthquake  loadings.  While  the  Westergaard  equation  for  approxi- 
mating the  hydrodynamic  interaction  between  the  dam  and  the  reservoir 
was  a major  contribution  in  its  time,  recent  studies  have  shown  that  ne- 
glecting the  compressibility  of  the  water  in  the  reservoir  can  lead  to 
serious  errors  in  estimating  the  natural  period  of  vibration  of  the  dam. 
The  effect  of  the  reservoir  water  is  to  cause  a significant  increase  in 
the  natural  period  of  vibration  of  the  dam.  This  lengthening  of  the 
natural  period  depends  on  the  depth  of  the  reservoir  water  and  the  modu- 
lus of  elasticity  of  the  dam.  For  a practical  range  of  values  of  elas- 
ticity and  a full  reservoir,  the  fundamental  period  may  increase  by  25 
to  50  percent. 6 


P.  Chakrabarti  and  A.  K.  Chopra,  "Earthquake  Analysis  of  Gravity  Dams 
Including  Hydrodynamic  Interaction,"  Earthquake  Engineering  and 
Structural  Dynamics,  Vol  2 (1973),  pp’ 143-lob. 


In  addition,  it  has  been  shown  that  hydrodynamic  forces  resulting 
from  the  vertical  component  of  earthquake  ground  motions  are  comparable 
to  those  caused  by  horizontal  earthquake  ground  motions  for  reservoirs 
of  moderate  depth,  but  much  larger  for  small  depths. 

Thus,  the  earthquake  factor  is  also  subject  to  significant  bias  and 
uncertainty. 


Bias  and  Uncertainty  in  the  Horizontal  Force  Model 

With  the  exception  of  wave,  ice,  wind,  and  horizontal  silt  pressure 
loading,  the  horizontal  force  model  provides  a reasonable  representation 
of  the  true  horizontal  forces  acting  on  the  dam  when  the  appropriate 
corrective  factors  are  introduced  into  the  equations  and  earthquake 
loadings  are  neglected.  The  bias  introduced  by  neglecting  the  wave, 
ice,  wind,  and  horizontal  silt  pressure  would  be  unconservative;  there- 
fore, a corrective  factor  greater  than  1 should  be  included  to  offset 
this  systematic  bias.  The  bias  and  uncertainty  can  only  be  assigned 
after  proper  consideration  is  given  to  the  magnitude  of  the  error  and 
its  potential  variability. 


Sliding  Resistance  Model 

Equations  for  Mean  and  Uncertainty 

The  sliding  resistance.  Q,  is  a function  of  the  unit  shear  strength 
at  no  normal  load  (cohesion),  the  angle  of  internal  friction  on  the  po- 
tential failure  surface,  and  the  vertical  forces  acting  on  the  dam,  in- 
cluding uplift,  as  expressed  in  Eq  3.  Thus,  in  computing  the  sliding 
resistance,  the  major  variables  are  s,  <J>,  H],  Hg,  and  K;  thus  corrective 
factors  Ns,  Na,  nH-|,  Nh2,  N|<,  respectively,  and  Ng  must  be  introduced 
into  Eq  3 to  Obtain  the  mean  value  of  the  total  sliding  resistance 
force,  Q.  This  substitution  yields 


LEq  26] 


The  total  uncertainty  associated  with  Q can  be  evaluated  through  Eq 
21.  Again,  if  only  Hi  and  H2  are  assumed  to  be  correlated  variables, 
the  total  uncertainty  becomes 
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Q * 'q  * Qf  |[Ytta,-KL)tan»]2  * [(y^W-y^ 

* l/2>BLK)tan»]2  N ^H2'^  * L2N2S2n2  * [l/2Tu(HrH2)]2 

’ "k'^K  * (Wl'VV  * ,/2-<„l1H3  • V2YwH2tan9"  [Eq  27] 

2 

- ywL[H2+l/2K(HrH2)sec24»j  N202Q2  + ^(L^KLjtan*] 

[(YwH2tan0"  - ywL  + l/2YwLK)tan4>]  P^h^h^h^I H2aH^H2 


Biur,  and  Uncertainty  of  the  Principal  Variables 

Shear  resistance  within  a foundation  and  between  a dam  and  its 
foundation  depends  on  the  cohesion  and  the  angle  of  internal  friction  in 
the  foundation  and  the  bond  between  the  concrete  and  foundation  at  the 
base  of  the  dam.  Traditionally,  the  shear  resistance  has  been  expressed 
as  varying  linearly  with  respect  to  the  normal  load  according  to  the 
Coulomb  equation,  the  variation  of  shear  resistance  versus  normal  load 
being  determined  by  laboratory  tests  and/or  in-situ  tests.  The  in-situ 
rock,  however,  will  invariably  contain  joints  or  other  types  of  geologi- 
cal discontinuities  which  may  be  spaced  from  just  a few  inches  to 
several  feet  apart.  Moreover,  a dam  will  usually  stress  a mass  of  rock 
large  enough  to  include  several  of  these  discontinuities.  Laboratory 
tests,  on  the  other  hand,  are  performed  on  intact  rock  cores  containing 
no  discontinuities  (or  at  most  a single  joint).  In  addition,  the  rock 
cores  used  in  the  laboratory  have  be  i removed  from  their  in-situ  envi- 
ronment, i.e.,  the  moisture  conditions  may  change,  the  confining  pres- 
sure has  been  altered,  and  the  procedures  of  sampling  may  have  degraded 
the  quality  of  the  sample.  A scale  effect  may  be  introduced  because  of 
the  size  of  the  sample  tested,  and  the  results  may  also  be  sensitive  to 
the  rate  of  testing. 

In-situ  testing  may  provide  a better  picture  of  the  quality  of  the 
rock,  but  it  is  still  indicative  rather  than  precise.  Each  test  pro- 
vides information  relative  to  one  small  area;  and,  in  general,  the  tests 
are  done  quickly  so  they  are  not  representative  of  the  behavior  of  the 
foundation  under  long  term  conditions. 
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Shear  Strength  (Cohesion).  The  pattern  of  joints,  shear  zones,  and 
faults  in  a rock  mass  make  the  shear  strength  of  rock  highly  aniso- 
tropic. When  the  direction  of  loading  is  such  that  the  particular  fail- 
ure surfaces  must  cut  across  the  structural  features,  the  shear  strength 
will  approach  that  of  the  intact  rock.  However,  when  the  direction  of 
loading  is  parallel  or  subparallel  to  the  structural  features,  the  shear 
strength  is  governed  by  the  shearing  resistance  along  the  rock  surfaces 
of  the  discontinuity;  and  the  effective  shear  strength  will  be  reduced 
to  a value  much  below  the  intact  rock  strength.  This  latter  case  is 
generally  recognized  as  being  critical  in  evaluating  the  sliding  re- 
sistance of  gravity  dams. 

The  shearing  resistance  of  rocks  has  been  shown  to  exhibit  both  a 
peak  shear  strength  and  a residual  strength  as  illustrated  in  Figure 
10.  ^ in  the  shearing  process,  the  shearing  strength  along  the  discon- 
tinuity reaches  a maximum  at  some  small  value  of  displacement  where 
fracture  occurs.  Thereafter,  the  shear  strength  gradually  decreases 
with  continuing  displacements  until  finally  at  large  displacement,  the 
shear  strength  approachs  asymptotically  the  residual  strength.  When 
more  failure  envelopes  are  drawn  through  the  maximum  and  minimum  values 
of  shear  strength  obtained  for  specimens  under  different  normal  loads, 

N,  the  results  may  be  similar  to  those  shown  in  Figure  11.  The  vertical 
distance  between  the  two  envelopes  indicates  the  reduction  in  shear 
strength  by  continued  displacement.  Note,  also,  that  the  residual 
strength  envelope  does  not  show  a cohesion  intercept  and  can  be  de- 
scribed uniquely  by  the  angle  of  residual  sliding  reaction.  Thus,  the 
shearing  resistance  along  a discontinuity  in  the  field  at  a given  value 
of  normal  load  is  dependent  on  the  magnitude  of  the  prior  relative  dis- 
placements which  have  occurred  between  the  rock  surfaces.  When  cohesion 
is  taken  into  account  in  the  design,  it  is  customary  to  adopt  values  1/3 
to  1/5  of  the  apparent  cohesion  indicated  by  representative  samples. 8 

Angle  of  Internal  Friction.  Typical  peak  and  residual  angles  of 
internal  friction  for  several  rock  types  are  presented  in  Table  3.  As  a 
result  of  this  variation  in  the  peak  and  residual  angle  of  internal 
friction,  the  coefficient  of  friction  (tan<j>)  along  the  failure  plane  can 
vary  by  a factor  ranging  between  1 and  2. 

Uplift.  A dam  is  subjected  to  water  pressure,  not  only  on  exposed 
faces  Dut  also  on  its  base.  The  pressure  on  the  base  (uplift)  reduces 
the  effective  weight  of  the  dam  and,  consequently,  affects  the  sliding 
resistance  of  the  dam.  The  intensity  of  the  hydrostatic  pressure  and 
the  area  upon  which  the  pressure  acts  have  been  the  subject  of 


7 K.  G.  Stagg  and  0.  C.  Zienkiewicz,  Rock  Mechanics  in  Engineering  Prac- 
o tice  (John  Wiley,  19681. 

^"Advances  in  Rock  Mechanics,"  Proceedings  of  the  Third  Congress  of  the 
International  Society  for  Rock  Mechanics,  Vol  I,  Part  A (National 
Academy  of  Sciences,  1^77). 
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Figure  10. 


Shear  strength  versus  horizontal  displacement  showing 
maximum  and  residual  stresses.  (Reprinted  with  per- 
mission from  K.  G.  Stagg  and  0.  C.  Zienkiewicz, 

Rock  Mechanics  in  Engineering  Practice  [John  Wiley, 


Figure  11.  Maximum  strength  and  residual  failure  envelope  for 
initially  intact  specimens.  (Reprinted  with  per- 
mission from  K.  G.  Stagg  and  0.  C.  Zienkiewicz, 

Rock  Mechanics  in  Engineering  Practice  [John  Wiley, 
1 968 JT 
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discussion  in  the  past;  however,  today  it  is  accepted  design  practice  to 
assume  the  pressure  acts  on  the  full  area  of  the  base. 


It  is  also  accepted  design  practice  to  adopt  a distribution  of 
uplift  pressure  as  shown  in  Figure  2,  the  value  of  K being  selected 
after  thorough  consideration  of  the  porosity  of  the  foundation  rock,  the 
existence  of  joints  and  cracks  therein,  and  the  effectiveness  of  the 
drains.  Based  on  experience,  the  U.S.  Bureau  of  Reclamation  normally 
uses  a value  of  K = 0.33;  the  Tennessee  Valley  Authority  has  used  0.25; 
while  the  Corps  of  Engineers  uses  values  between  0.25  and  0.50, 
depending  on  what  is  considered  appropriate  for  the  site. 9 

In  summary,  significant  bias  and  uncertainty  is  associated  with  se- 
lection of  the  cohesion,  angle  of  internal  friction,  and  uplift.  In  ad- 
dition, additional  effects  that  were  not  discussed  may  increase  this 
bias  and  uncertainty,  such  as  rate  of  testing  and  sample  size. 


Bicu;  and  UncnrbainLij  in  tin'  dlidin.j  iu-nic-tarw  M>./W 

The  relationship  for  the  sliding  resistance  neglects  the  effects  of 
vertical  earthquake  acceleration,  which  is  typically  taken  as  two-thirds 
of  the  horizontal  earthquake  acceleration  for  design  purposes.  The 
gross  effect  of  the  vertical  acceleration  is  to  reduce  the  effective 
normal  force  acting  on  the  interface  at  the  base  of  the  dam  and  thus 
reduce  the  sliding  resistance. 


Table  3 

Maximum  and  Residual  Angles  of  Internal  Friction 
(From  K.  G.  Stagg  and  0.  C.  Zienkiewicz,  Rock  Mechanics  in  Engineering 
Practice  (John  Wiley,  1968). 


Maximum 

Residual 

Rock  Type 

? 

* 

(deg) 

(deg) 

Sandstone 

48-50 

25-34 

Limestone 

37-58 

33-37 

Quartz  Monzonite 

56 

32 

9 H.  H.  Thomas,  The  Engineering  of  Large  Dams,  Part  I (John  Wiley, 
1976),  p 257. 
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Moreover,  since  the  shear  stress  in  the  foundation  of  dams  is  usu- 
ally not  distributed  uniformly  along  the  potential  failure  surfaces, 
progressive  failure  may  occur.  This  failure  consists  of  mobilization  of 
the  peak  shear  strength  in  certain  zones  of  the  failure  surface,  while 
other  zones  have  fractured  and  the  shear  strength  is  gradually 
decreasing.  Within  the  framework  of  present  knowledge,  it  is  not  possi- 
ble to  accurately  calculate  the  maximum  shearing  resistance  that  may  be 
developed  during  progressive  collapse;  however,  it  is  possible  to  bound 
the  problem.  An  upper  limit  corresponds  to  the  mobilization  of  the  peak 
shear  strength  of  the  material  over  the  potential  failure  surface,  and 
the  lower  limit  corresponds  to  the  mobilization  of  the  residual  shear 
strength.  For  reasons  of  safety,  the  lower  limit  is  usually  taken  as 
the  strength  along  the  potential  rupture  surface. 


Foundation  Pressure  Model 


Equations  for  Mean  and  Uncertainty 


The  maximum  foundation  pressure,  including  uplift,  may  be  obtained 
from  Eq  14.  By  inspection,  the  foundation  pressure  is  a function  of  the 
variables  , H2,  a,  C,  and  K.  Introducing  the  appropriate  corrective 
factors  and  dropping  the  subscripts  and  superscripts  associated  with  P- 

*ields  r-r  o o 

A,L<VrH3> + '/Vj  + '''2NH2Wa"e" 


P = N, 


{{v 


-VtNH2V'/2NKK<HH)HrNH2H2)]}i’/u 

* {[4/,5Wh,“CH1  * ,/6\.NH,H?  - ,/6>„NH2H2  Uq  26] 


* Naavc(l/6L,H2+l/2L2H2tl/6L3H2)  - ) 


- 1/2ywL1H3(1/2L-1/3L1  ) - l/2ycL1H3(l/2L-2/3L1)  - Yc4Ho(l/2L 
-L1-1/2L2)  + 1/2ycL3H4(1/2L-2/3L3)  + l/2YcN^tanO' 

• ( 1 /2L- 1 / 3NH^H2tan0 ’ ) + 1/12ywNkK(Nh^H1-Nh^H2)L2](1/L2) 

+ TwNH2H2f  secte"  " TwNh  H2^sin0"+COsO"^j 


The  total  uncertainty  associated  with  P can  be  evaluated  through  Eq 
21.  To  improve  the  tractability  of  the  derivation,  the  expressions  for 
the  derivatives,  evaluated  at  their  mean  values,  are  included.  Again, 
the  subscripts  and  superscripts  associated  with  P have  been  dropped. 

Ufj-"  jyw(Lrl/2KL)(l/L)  + 6[12/15CctH^+l/2YwH^-l/2ywL1  (L-L] ) 


-1/6y, 


wKL2](1/L2)] 


sec2e"  = B, 


= |rw(‘l-i+H2tan0"-L+l/2K)(l/L)  + 6yw[-l/2H2+l/2LH2tane" 

-l/2H2tano"  - 1/6KL2](1/L2)  + Yl~|sec2e"  - yJsinQ" 

^ W I w 


+ cos0")tan2e"  = B. 


[Eq  29] 


= 6(4/15ctH3)(l/L2)  = B3 


SK  = •3/2yw(H1'H2)  = B4 


1 5 = 6[4/15CHj  + Yc(1/6L1H^1/2L2H^1/6L3H2)](1/L2)  = B& 


Therefore,  the  total  uncertainty  of  the  maximum  foundation  pressure 
becomes 


% = V ^ + B2NH2H2«H2  + B3NcC^ 

+ b4nkk^k  + 


[Eq  30] 


1 * n -f  1 1 ^iT1' - ' - - 


Bias  and  Uncertainty  of  the  Principal  Variables 

The  foundation  pressure  is  a function  of  the  random  variables  pre- 
viously discussed;  consequently,  the  discussions  regarding  the  bias  and 
uncertainty  of  each  of  these  random  variables  is  also  applicable  to  the 
foundation  pressure. 


Bias  and  Uncertainty  of  the  Foundation  Pressure  Model 


In  the  gravity  method  of  analysis,  each  section  of  the  dam  is  as- 
sumed to  act  independently,  and  the  loads  are  transmitted  to  the  founda- 
tion by  cantilever  action.  In  this  method,  the  effects  of  elasticity  of 
the  dam,  as  well  as  the  foundation,  are  ignored.  In  reality,  both  the 
dam  and  its  foundation  are  elastic,  and  the  effects  of  this  elasticity 
will  cause  the  foundation  pressure  to  depart  from  the  linear  distribu- 
tion inherently  assumed  through  the  use  of  Eqs  11  and  12.^0 

For  example.  Figure  12  illustrates  the  foundation  pressures  mea- 
sured under  the  Shasta  Dam,  in  January  1954  at  the  time  the  -eservoir 
was  nearly  full.  The  measured  maximum  stress  was  4.1  MPa  at  a point 
near  the  center  of  the  base.  Near  the  downstream  face,  the  foundation 
pressure  reduces  to  about  1 MPa,  which  indicates  the  downstream  third  of 
the  concrete  dam  was  carrying  a relatively  minor  part  of  the  load.  The 
maximum  measured  shearing  stress  occurred  near  the  downstream  face  in 
the  area  of  minimum  vertical  pressure.  Measurements  at  Grand  Coulee  Dam 
show  a quite  different  pattern  (Figure  13).  Thus,  there  may  be  consid- 
erable bias  and  uncertainty  in  the  foundation  pressures  model. 


Foundation  Compressive  Strength  Model 

Equations  for  Mean  and  Uncertainty 

Usually  the  compressive  strength  of  the  foundation  material,  o, 
will  be  determined  from  laboratory  tests  on  small  samples.  However, 
mature  judgment  must  supplement  the  laboratory  test  results  in  order  to 
estimate  the  compressive  strength  of  the  in-situ  foundation  material, 
i.e.,  evaluate  the  effects  of  lateral  pressures,  cracks,  and  joints, 
solubility,  etc.  Thus,  to  estimate  the  mean  compressive  strength  of  the 
foundation  material,  corrective  factors  Na  and  Nu  are  introduced  such 
that  the  mean  value  of  the  compressive  strength  becomes 


U = N (N 
u'  a 


o) 


[Eq  31] 


10  H.  H.  Thomas,  The  Engineering  of  Large  Dams,  Part  I (John  Wiley, 
1976),  p 257. 
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Figure  12. 


Foundation  pressure  measurements  under  Shasta  Dam. 
(Reprinted  with  permission  from  H.  H.  Thomas, 
TfteEngj peering  of  Large  Dams.  Part  I (John  Wiley, 
1976; , p 257. 


Figure  13.  Foundation  pressure  measurements  under  Grand  Coulee  Dam 

(Reprinted  with  permission  from  H.  H.  Thomas,  The  Enaineer 
T_n9  of  Large  Dams,  Part  I (John  Wiley,  1976),  p~257. 


Since  the  mean  value  is  independent  of  the  other  major  variables, 
the  total  uncertainty  associated  with  it  becomes 


n?  = 5/  + y [N202ft>2] 

u u jjZ  L a a 


[Eq  32] 


Unomrtainty  in  Foundation  Contprea&im  strength 


The  strength  of  the  foundations  on  which  dams  are  founded  varies 
from  solid  rock,  which  is  much  stronger  than  concrete,  through  all 
stages  of  decomposed  and  disintegrated  rock.  Yet,  the  compressive 
strength  of  the  Toundation  rock  is  very  important  in  determining  the 
base  width  of  the  dam.  Normally,  the  maximum  allowable  compressive 
stress  in  the  foundations  is  25  percent  and  80  percent  of  the  foundation 
compressive. strength  for  normal  and  extreme  loading  conditions,  re- 
spectively. ' 


Bias  and  Uncertainty  of  Comgrcsaive  Strength  Model 

If  a proper  corrective  factor  is  applied  to  the  foundation  com- 
pressive strength,  then  the  bias  associated  with  the  model  is  unity. 
The  uncertainty  is  the  confidence  one  has  that  all  corrective  factors 
are  accounted  for  properly. 


11 


A.  R.  Golze,  Handbook  of  Dam  Engineering  (Van  Nostrand  Reinhold  Com- 
pany, 1977). 


5 SIMPLE  APPLICATIONS 


1 


i 


I 


General 


In  the  interest  of  simplicity  and  for  the  lack  of  design  infor- 
mation and  data  for  an  actual  dam,  examples  of  two  moderately  low  grav- 
ity dams  designed  in  accordance  with  conventional  design  procedures  were 
selected  for  analysis  to  illustrate  the  probabilistic  concepts  developed 
in  the  preceding  section.  These  example  designs,  which  are  contained  in 
Engineering  for  Dams,  were  designed  to  have  a minimum  shear-friction 
safety  factor  of  5.0  and  an  overturning  safety  factor  of  1.0. ^ in  one 
example,  a 200-ft  nonoverflow  gravity  dam  was  designed  considering  only 
water  loadings;  in  the  other  example,  an  identical  dam  was  designed  con- 
sidering both  water  and  earthquake  loadings.  The  primary  objectives  of 
the  analyses  in  this  chapter  are  to  estimate  the  probability  of  failure 
of  the  dams  with  respect  to  sliding  and  overturning  and  to  develop  an 
insight  into  the  effects  of  various  sources  of  uncertainty  underlying 
the  design  parameters  on  the  safety  of  the  dam. 


200-Ft  Nonoverflow  Gravity  Dam 

The  cross  section  of  the  200-ft  nonoverflow  gravity  dam,  which  was 
designed  to  resist  only  water  loadings,  is  shown  in  Figure  14  together 
with  the  design  specifications  and  conditions. 

Since  the  statistics  of  the  principal  random  variables  were  not 
given  in  the  cited  example,  it  was  necessary  to  estimate  reasonable 
values  for  the  lower  and  upper  bounds  of  the  mean,  to  assume  some  appro- 
priate distributions  for  the  mean  over  its  range,  and  to  compute  the 
mean  and  coefficient  of  variation  by  standard  statistical  methods.  For- 
mulas for  computing  the  mean  and  coefficient  of  variation  for  some 
simple  distributions  are  presented  in  Table  4.^3  Next,  estimates  for  the 
corrective  factor  associated  with  each  variable  were  developed,  based  on 
the  considerations  discussed  in  the  preceding  section.  A reasonable 
value  was  assigned  to  each  of  the  corrective  factors,  then  a series  of 
lower  and  upper  values  was  established  in  an  effort  to  bound  the  likely 
corrective  factors.  The  coefficient  of  variation  for  the  corrective 
factors  was  assigned  a value  of  0.2,  based  on  intuition  and  professional 
judgment.  Finally,  estimates  of  the  bias  associated  with  the  horizontal 
force,  sliding  resistance,  foundation  pressure,  and  foundation  com- 
pressive strength  models  were  developed.  A mean  value  of  1.05  was  as- 


12W.  P.  Creager,  J.  D.  Justin,  and  J.  Hinds,  Engineering  for  Dams  (John 
Wiley,  1945). 

13M.  S.  Yuceman,  W.  H.  Tank,  and  A.  H-S.  Ang,  A Probabilistic  Study  of 
Safety  and  Design  of  Earth  Slopes,  Civil  Engineering  Studies  Structural 
Research  series  No.  402  (university  of  Illinois,  1973). 
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Maximum  Depth  of  Water  to  Be  Retained  = 200  ft 
Depth  of  Tailwater  = 0 
Crest  Width  = 24  ft 

Spill  Way  Crest  to  Maximum  Water  Surface  = 10 
Weight  of  Concrete  = 150  pcf 

Angle  of  Internal  Friction  for  Foundation  Material  = 
Ultimate  Shear  Resistance  of  Foundation  = 800  psi 
Maximum  Allowable  Inclined  Stress  = 60  ksf 
Weight  of  Water  =62.5  pcf 
Uplift  Intensity  Factor  =0.5 

Figure  14.  200-ft  nonoverflow  gravity  dam 


Mean  and  Coefficient  of  Variation  of  Nj  Corresponding  to  Different  Distributions  Assumed  Over  Its  Range. 

From  M.  S.  Yuceman,  W.  H.  Tank,  and  A.  H-S.  Ang,  A Probabilistic  Study  of  Safety  and  Design  of  Earth  Slopes, 
Civil  Engineering  Studies  Structural  Research  Series  No.  402  (University  of  Illinois,  1973). 


lower  bound  of  N 


signed  to  the  horizontal  force  because  factors  such  as  wind,  wave,  ice, 
and  silt  pressure  were  neglected  in  computing  the  horizontal  force 
acting  on  the  dam.  A mean  value  of  0.90  was  assigned  to  the  sliding  re- 
sistance in  partial  recognition  of  potential  progressive  failure.  In 
the  absence  of  more  definitive  data,  a mean  value  of  1.0  was  assigned  to 
both  the  foundation  pressure  and  compressive  strength.  The  coefficient 
of  variation  for  all  the  bias  factors  was  assigned  a value  of  0.2.  The 
mean  values  and  coefficients  of  variations  for  the  principal  variables, 
corrective  factors,  and  bias  are  summarized  in  Table  5. 

In  the  analysis,  the  mean  value  of  each  corrective  factor  was  sys- 
tematically varied,  and  the  probabilities  of  sliding  and  overturning 
failures  were  computed.  For  example,  in  computing  the  probability  of 
sliding  failure,  the  corrective  factor  for  each  of  the  four  variables 
influencing  a sliding  failure  was  systematically  varied  while  the  cor- 
rective factors  for  the  other  variables  were  assigned  their  intermediate 
mean  value.  That  is,  when  the  corrective  factor  for  a variable  such  as 
the  reservoir  water  elevation  was  varied,  the  corrective  factors  for  the 
angle  of  internal  friction,  hydrostatic  pressure  intensity  factor,  and 
foundation  shear  strength  were  1.5,  1.0,  and  2.0,  respectively.  In  this 
fashion,  the  effect  of  a particular  parameter  could  be  isolated  to  the 
maximum  extent.  A similar  procedure  was  followed  in  computing  the 
probability  of  an  overturning  failure. 

The  probability  of  sliding  failure  of  the  200-ft  nonoverflow  dam 
was  less  than  10“6  for  all  cases  and  outside  the  range  of  practical  in- 
terest. Thus,  detailed  results  are  not  included.  On  the  other  hand, 
the  probability  of  overturning  failure  was  greater  than  10”^  as  illus- 
trated in  the  plots  of  probability  of  overturning  failure  versus  the 
corrective  factors  for  reservoir  water  elevation,  hydrostatic  pressure, 
intensity,  and  foundation  compressive  strength  shown  in  Figures  15,  16, 
and  17,  respectively.  While  the  probability  of  overturning  failure 
might  initially  appear  to  be  insensitive  to  the  corrective  factor  for 
reservoir  water  elevation,  it  should  be  noted  that  the  corrective  factor 
is  only  varied  through  a range  of  plus  and  minus  4 percent,  and  this 
small  change  results  in  nearly  an  order  of  magnitude  difference  in  the 
computed  probability  of  overturning  failure.  The  corrective  factor  for 
hydrostatic  pressure  intensity  has  virtually  no  effect  on  the  proba- 
bility of  overturning  failure  for  the  range  considered.  As  one  would 
expect,  the  corrective  factor  for  foundation  compressive  strength  has  a 
great  effect  on  the  probability  of  overturning  failure.  Increasing  the 
corrective  factor  from  1.0  to  4.0  decreases  the  probability  of  over- 
turning failure  from  nearly  10-1  to  10‘5,  i.e.,  about  four  orders  of 
magnitude.  Thus,  the  corrective  factor  for  foundation  compressive 
strength  should  be  established  with  care  in  any  analysis  of  an  actual 
dam. 
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Figure  15.  Probability  of  overturning  failure  vs. 

reservoir  water  elevation  corrective  factor. 
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HyOroatatic  Preaeure  In  family  Corrective  factor 


Figure  16.  Probability  of  overturning  failure  vs. 

hydrostatic  pressure  intensity  corrective  factor. 
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Figure  17.  Probability  of  overturning  failure  vs.  foundation 
compressive  strength  corrective  factor. 
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200-Ft  Nonoverfl ow  Gra v j_ly  _0am_ W i th  Earthquake  Loading 

The  cross  section  of  the  200-ft  nonoverflow  gravity  dam  designed  to 
resist  a moderate  to  severe  earthquake  in  addition  to  water  loading  is 
shown  in  Figure  13,  together  with  the  design  specifications  and  condi- 
tions. 

Again,  since  the  statistics  of  the  principal  random  variables  were 
not  presented  in  the  cited  example,  it  was  necessary  to  follow  a pro- 
cedure similar  to  that  employed  in  the  previous  application  to  develop 
means  and  coefficients  of  variation.  However,  statistics  for  the 
earthquake  acceleration  coefficient  and  the  earthquake  factor  had  to  be 
included.  These  were  assumed  to  have  mean  values  of  0.1  and  51.7,  re- 
spectively, and  coefficients  of  variations  0.3  and  0.005,  respectively. 

A summary  of  the  means  and  coefficients  of  variation  for  the  principal 
variables,  corrective  factors,  and  bias  parameters  is  presented  in  Table 
6. 

For  the  200-ft  nonoverflow  dam  with  both  water  and  earthquake  load- 
ings, the  probability  of  sliding  failure  was  less  than  10"'*  for  all 
cases.  While  these  probabilities  are  outside  the  range  of  interest,  an 
indication  of  how  the  computed  probability  of  sliding  failure  varies 
with  the  earthquake  acceleration  coefficient  and  foundation  shear 
strength  corrective  factors  can  be  derived  from  Figures  19  and  20.  With 
the  exception  of  these  corrective  factors,  the  probability  of  sliding 
failure  was  relatively  insensitive  to  variations  in  the  other  corrective 
factors  and,  in  general,  the  probability  of  sliding  failure  was  in  the 
order  of  10  . 


^ The  probability  of  overturning  failure,  however,  was  greater  than 
10"  for  all  cases  and,  in  general,  was  in  the  order  of  10”^.  The 
effect  of  the  various  corrective  factors  on  the  computed  probability  of 
overturning  failure  can  be  observed  in  Figures  21  through  25. 

Variations  in  the  foundation  compressive  strength  and  earthquake  accel- 
eration coefficient  corrective  factors  produced  the  greatest  impact  on 
the  computed  probability  of  overturning  failure.  Thus,  special  care 
should  be  observed  in  assigning  values  for  these  corrective  factors. 


The  results  cited  so  far  primarily  illustrate  the  effects  that 
variations  in  the  corrective  factor  have  on  the  probability  of  a sliding 
or  overturning  failure,  i.e.,  the  effects  of  varying  the  corrective 
factor  uncertainty,  bias,  and  bias  uncertainty  have  not  been  in- 
vestigated. To  evaluate  the  effect  that  corrective  factor  uncertainty 
may  have  on  the  probability  of  failure,  the  analysis  required  to  obtain 
the  data  presented  in  Figure  25  was  repeated  for  a series  of  cases  in 
which  the  corrective  factor  uncertainty  was  assigned  values  of  0,  0.1, 
0.2,  and  0.3.  These  results  are  summarized  in  Figure  26,  and  they  indi- 
cate that  corrective  factor  uncertainty  becomes  more  significant  as  the 
magnitude  of  the  corrective  factor  increases.  Moreover,  for  the  most 
likely  range  of  foundation  compressive  strength  corrective  factor  (2-4), 
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Maximum  Depth  of  Water  to  Be  Retained  = 200  ft 
Depth  of  Tailwater  = 0 
Crest  Width  = 24  ft 

Spill  Way  Crest  to  Maximum  Water  Surface  = 10  ft 
Uplift  Intensity  Factor  =0.5 

Angle  of  Internal  Friction  for  Foundation  Material  = 36.9° 

Ultimate  Shear  Resistance  of  Foundation  = 115.2  ksf 

Maximum  Allowable  Inclined  Stress  = 60  ksf 

Ratio  of  Earthquake  Acceleration  to  Gravity  =0.1 

Earthquake  Factor  = 51.7 

Weight  of  Concrete  = 150  pcf 

Weight  of  Water  =62.5  pcf 

Figure  18.  200-ft  nonoverflow  gravity  dam  with  earthquake  loading. 
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. Probability  of  sliding  failure  vs.  earthquake 
acceleration  coefficient  corrective  factor. 


Figure  20.  Probability  of  sliding  failure  vs.  foundation 
shear  strength  corrective  factor. 


Reeorvok  Water  Elevation  Corrective  Factor 


servolr  water  elevation 


Figure  21.  Probability  of  overturning  failure  vs 
corrective  factor. 


Figure  22.  Probability  of  overturning  failure  vs.  earthquake  acceleration 
coefficient  corrective  factor. 


Earthquake  Factor  Corrective  Factor 

Figure  23.  Probability  of  overturning  failure  vs.  earthquake  factor 
corrective  factor. 


Hydrostotic  Pressure  Intensity  Corrective  Factor 

Figure  24.  Probability  of  overturning  failure  vs.  hydrostatic  pressure 
intensity  corrective  factor. 


fbaadotiaw  CMimiH  Slrtngtk  CorracH**  Factor 

Probability  of  overturning  failure  vs.  foundation  compressive 
strength  corrective  factor. 


Figure  25 


Figure  26.  Effect  of  corrective  factor  uncertainty 


the  difference  between  assigning  a value  of  0 or  0.3  to  the  corrective 
factor  uncertainty  will  change  the  results  by  about  an  order  of  mag- 
nitude. 

Likewise,  to  evaluate  the  effect  that  variations  in  the  bias  may 
have  on  the  probability  of  failure,  a similar  analysis  was  repeated  for 
a series  of  cases  in  which  the  bias  for  the  foundation  pressure  and  the 
foundation  compressive  strength  were  jointly  assigned  values  of  0.8, 

1.0,  1.2,  and  1.4.  While  the  mean  values  for  the  load  and  resistance 
varied  in  accordance  with  the  assigned  values  of  bias,  the  computed 
probability  of  overturning  did  not  depart  from  the  results  presented  in 
Figure  25  in  sufficient  magnitude  to  warrant  plotting.  Thus,  establish- 
ing the  proper  ratio  between  the  bias  for  the  foundation  pressure  and 
the  foundation  compressive  strength  is  probably  more  important  than 
knowing  the  precise  magnitude  of  the  bias. 

Finally,  a similar  analysis  was  repeated  to  evaluate  the  effect  of 
bias  uncertainty.  Values  of  0,  0.1,  0.2,  and  0.3  were  assigned  to  the 
bias  uncertainty,  and  as  before,  the  probability  of  overturning  failure 
was  computed.  These  results  are  presented  in  Figure  27  and,  in  general, 
they  are  similar  to  the  results  presented  in  Figure  26.  However,  these 
results  indicate  that  the  bias  uncertainty  has  a slightly  greater  effect 
on  the  probability  of  failure  than  the  corrective  factor  uncertainty. 
Thus,  as  larger  values  are  assigned  to  the  corrective  factors,  greater 
care  should  be  observed  first  in  assigning  values  to  the  bias  uncer- 
tainty and  second  in  assigning  values  to  the  corrective  factor  uncer- 
tainty. 


6 PRELIMINARY  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


Concepts  and  methods  of  probability,  i.e.,  extended  reliability 
theory,  provide  the  proper  basis  for  evaluating  of  the  safety  of  dams. 
Moreover,  implementation  of  these  concepts  and  methods  is  within  the 
current  state  of  the  art  as  illustrated  by  the  probabilistic  concept  for 
concrete  gravity  dams  summarized  herein.  With  such  formulations,  the 
effects  of  different  sources  of  uncertainty  underlying  the  conventional 
design  parameters  can  be  systematically  combined  and  analyzed  in  a 
manner  suitable  for  a quantitative  assessment  of  the  safety  of  a dam 
with  respect  to  various  modes  of  failure.  Such  formulations  also  pro- 
vide an  improved  basis  for  updating  design  criteria.  Finally,  it  should 
be  emphasized  that  the  lack  of  statistical  data  is  not  a valid  reason 
for  rejecting  probabilistic  concepts  and  methods.  Indeed,  it  is  only 
through  probabil istic  models  that  the  significance  of  objective  infor- 
mation, or  lack  thereof,  can  be  properly  assessed  and  combined  with 
engineering  judgment  to  provide  a quantitative  assessment  of  the  safety 
of  a dam. 

For  the  simple  applications  of  the  probabilistic  concept  for  con- 
crete gravity  dam  analysis  presented  herein,  it  was  observed  that: 

1.  The  probability  of  sliding  failure  was  most  sensitive  to  the 
corrective  factors  assigned  to  the  foundation  shear  strength  and  the 
earthquake  acceleration  coefficient.  There  was  considerably  less  sensi- 
tivity associated  with  the  other  corrective  factors. 

2.  The  probability  of  overturning  failure  was  most  sensitive  to 
the  corrective  factors  assigned  to  the  foundation  compressive  strength 
and  the  earthquake  acceleration  coefficient.  The  corrective  factors  as- 
signed to  reservoir  water  elevation  and  the  earthquake  factor  produced 
considerably  less  sensitivity  in  the  results,  and  the  results  were  vir- 
tually insensitive  to  the  hydrostatic  pressure  intensity  corrective 
factor.  The  probability  of  overturning  failure  also  became  more  sensi- 
tive as  large  values  were  assigned  to  the  corrective  factor  for  bias  un- 
certainty and  corrective  factor  uncertainty,  with  the  bias  uncertainty 
having  a slightly  greater  effect.  Moreover,  establishing  the  proper 
ratio  between  the  bias  associated  with  the  load  and  resistance  models 
may  be  more  important  than  determining  the  precise  magnitude  of  the 
bias. 


3.  In  any  future  application  of  this  probabilistic  concept  for 
gravity  dam  analysis,  special  care  should  be  exercised  in  assigning 
values  to  the  corrective  factors  for  the  foundation  shear  strength,  the 
foundation  compressive  strength,  the  earthquake  acceleration  coeffi- 
cient, the  ratio  of  the  bias  in  the  load  and  resistance  models,  and  — 
when  the  corrective  factors  are  large  — the  bias  uncertainty  and  cor- 
rective factor  uncertainty. 
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4.  For  the  design  example  where  only  water  loading  was  considered, 
the  probability  of  sliding  failure  was  less  than  10*6,  and  the  proba- 
bility of  overturning  failure  was  greater  than  10”5  for  the  range  of  pa- 
rameters considered. 

5.  For  the  design  example  where  both  water  and  earthquake  loadings 
were  considered,  the  probability  of  sliding  failure  was  less  than  10’^ 
and  the  probability  of  overturning  failure  was  greater  than  10~3  for  the 
range  of  parameters  considered. 


Recommendations 


The  results  to  date  illustrate  that  it  is  possible  to  quantify  the 
safety  of  a dam  in  a probabilistic  sense.  However,  the  results  and  con- 
clusions cited  above  should  be  viewed  as  preliminary  until  the  concept 
can  be  discussed  with  Corps  designers  experienced  in  gravity  dam  design, 
and  until  the  safety  of  a Corps  dam  can  be  analyzed  in  cooperation  with 
Corps  designers  to  provide  a practical  application  for  further  evalu- 
ation of  the  probabilistic  concept.  Finally,  the  concept  must  be  trans- 
lated into  a practical  format  that  can  be  readily  used  by  Corps  design 
engineers. 
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APPENDIX:  LIST  OF  SYMBOLS 


H1  ,H2 
H3,H4’ 


K 

L 


= area  of  base  of  the  dam  for  a section  of  unit  length 

= earthquake  factor  (Figure  3) 

= cumulative  distribution  function  of  n^ 

= summation  of  all  horizontal  forces  acting  on  dam  per 
unit  length 

= total  height  of  dam  section 

= height  of  reservoir  and  tailwater,  respectively 

= height  of  changes  in  section  in  the  upstream  and  down- 
stream faces  of  dam,  respectively 

= moment  at  inertia  of  a unit  length  of  the  base  about 
its  center  of  gravi ty 

= hydrostatic  pressure  intensity  factor 

= total  width  of  the  dam  at  the  base 


L-j  ,L2 »L3  = width  of  various  sections  of  the  dam 

M = overturning  moment,  including  uplift  moment,  about  the 

0 center  of  gravity  of  the  base 

m = distance  from  the  center  of  gravity  of  the  base  of  the 

dam 


N 


= correction  factor  with  mean  N and  coefficient  of  vari- 
ation A 


P-j  = total  reservoir  water  force  per  unit  length 

P^  = total  tailwater  force  per  unit  length 

P = total  hydrodynamic  force  due  to  earthquake  per  unit 

e length 

PQ  = total  inertial  force  due  to  earthquake  per  unit  length 

P.',P."  = inclined  pressure  at  the  upstream  and  downstream  faces, 

respectively 

P ' ,P  " = external  normal  unit  pressures  on  upstream  and  downstream 

force  of  the  dam,  respectively 


66 


p ' p " 

r ’ r 


P ' ,P  " 
v v 


P ' ,P  " 
u u 


R,S 

R,S 

R,S 

s 

Ss-f 

le 

U 

£V 

W 

o 

wrw2,w3 

Wu 

Xi 

a 

A 


Xi 


ar,as 


= vertical  pressure  at  the  upstream  and  downstream  faces, 
respectively 

= maximum  vertical  unit  pressures  on  upstream  and  down- 
stream faces  of  the  dam,  respectively 

= unit  uplift  pressures  on  the  base  of  the  dam  on  the 
upstream  or  downstream  force  of  the  dam,  respectively 

= probability  of  failure 

= sliding  resistance  of  dam  per  unit  length 

= true  random  resistance  and  true  random  load  (or  load 
effects),  respectively 

= models  of  R and  S 

A A 

= mean  values  of  R and  S 

= unit  shearing  st  ess  along  potential  failure  plane 

= shear  friction  safety  factor 

= period  of  vibration  of  the  earthquake 

= compressive  strength  of  the  foundation  material 

= summation  of  all  vertical  forces  acting  on  dam,  in- 
cluding uplift  but  excluding  foundation  reaction,  per 
unit  length 

= total  dead  weight  of  the  dam  per  unit  length 

= vertical  weight  of  reservoir  and  tail  water  component 
acting  on  dam  per  unit  length 

= total  uplift  force  per  unit  length 

= basic  resistance  or  load  variables 

= ratio  of  the  earthquake  acceleration  to  gravity 

N 

= coefficient  of  variation  of  representing  a measure 
of  the  errors  in  predicting  X. 

= coefficients  of  variation,  representing  a measure  of 
the  errors  in  predicting  R and  S,  respectively 


T 


@ 


■i 


* 
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x. 

Ai 


6r,6s 


°R’°S 


O'  ,0" 


nh,fiq,fip,Jo 


^R’^S 


= coefficient  of  variation  of  X.,  representing  a measure 
of  the  basis  variability  in  xj 


= coefficients  of  variation,  representing  a measure  of  the 
basic  variabilities  in  R and  S,  respectively 


= unit  compressive  strength  of  the  foundation  material 
= standard  deviation  of  R and  S,  respectively 


= angle  of  inclination  with  the  vertical  of  the  upstream 
and  downstream  faces  of  the  dam,  respectively 


= angle  of  internal  friction  for  the  foundation  material 


= coefficients  of  variation,  representing  the  measures  of 
the  total  uncertainties  in  the  theoretical  functions, 
respectively 


= coefficient  of  variation,  representing  measure  of  the 
total  uncertainties  in  X. 


'1 


coefficients  of  variation,  representing  measures  of  the 
total  uncertainties  in  R and  S,  respectively 


= correlation  coefficient 
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